Density distribution of a trapped two-dimensional strongly interacting Fermi gas 
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We calculate and measure the density distribution and cloud size of a trapped two-dimensional ®Li 
Fermi gas near a Feshbach resonance at low temperatures. Density distributions and cloud sizes are 
calculated for a wide range of interaction parameters using a local density approximation (LDA) and 
a zero-temperature equation of state obtained from quantum Monte Carlo simulations reported by 
G. Bertaina and S. Giorgini, Phys. Rev. Lett. 106, 110403 (2011). We find that LDA predictions 
agree well with experimental measurements across a Feshbach resonance. Theoretical results for 
Tan's contact parameter in a trapped gas are reported along with predictions for static structure 
factor at large momentum which could be measured in future Bragg spectroscopy experiments on 
two-dimensional Fermi gases. 
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I. INTRODUCTION 



The striking ability to manipulate and control ultra- 
cold atomic ^Li and ''"K Fermi gases has allowed the 
experimental investigation of strongly interacting two- 
component Fermi gases P, 0|. Using highly anisotropic 
pancake-shaped potentials to confine atoms in the lowest 
axial mode [3|-l5| , it has become possible to realize exper- 
imentally two-dimensional (2D) Fermi systems. At low 
temperatures these can exhibit exotic properties such as 
Berezinskii-Kosterlitz-Thouless (BKT) or inhonio- 

geneous Fulde-Ferrell-Larkin-Ovchinikov (FFLO) super- 
fluidity (9l4T^ . It has also been proposed that a 2D inter- 
acting Fermi gas may provide useful insights into high- 
temperature superconductivity (Tsl . [l^ and itinerant fer- 
romagnetism [isl. [l6j. To date, a weakly interacting 2D 
Fermi gas has been imaged in situ Q and used to char- 
acterize the crossover from two to three-dimensions Q. 
The observation of 2D confinement induced resonances 
and measurement of the molecular binding energy using 
rf-spectroscopy of a strongly interacting Fermi gas was 
recently reported However, the thermodynamic prop- 
erties of strongly interacting Fermi gases and fermionic 
mixtures [l7| in reduced dimensions are yet to be fully 
explored. 

On the theoretical side, numerous studies of 2D Fermi 
gases have been presented, addressing superfluid transi- 
tions 0, the effects of harmonic trap ping (l9| . and 
population and mass imbalance d, H^Hfj. Of particular 
importance, the equation of state of a 2D uniform Fermi 
gas was recently obtained through quantum Monte-Carlo 
(QMC) simulations through the crossover from a 
Bardeen-Cooper-SchriefFer (BCS) superfluid of Cooper 
pairs to a Bose-Einstein condensate (BEC) of tightly- 
bounded molecules Tan's universal many-body con- 
tact parameter was also found in 2D using the adiabatic 
relation [2^ . While the most theoretical studies have 
relied on a perturbative or mean-field approach, the ab- 



initio QMC results at zero temperature [23 should pro- 
vide a quantitative description of the many-body ground 
state of a strongly interacting 2D Fermi gas. 

In this work, we measure the density distribution and 
cloud size of a two-dimensional trapped ®Li Fermi gas at 
low temperatures in the BEC-BCS crossover and com- 
pare the data with theoretical predictions. The theoret- 
ical density distribution is calculated using a local den- 
sity approximation (LDA) based on the zero tem- 
perature QMC equation of state We find good 
qualitative agreement between experiment and theory in 
the strongly interacting regime near Feshbach resonance. 
When the atom number becomes large, we observe also 
substantial deviation from the 2D equation of state when 
new transverse vibrational modes are populated. 

We also give theoretical predictions for Tan's contact 
parameter in a trapped gas and the static structure fac- 
tor at large momentum, which could be measured in fu- 
ture Bragg spectroscopy experiments on 2D Fermi gases. 
Nontrivial pair correlations are reflected in the many- 
body part of the contact parameter. The many-body 
contribution to the contact exhibits a maximum near the 
3D Feshbach resonance. In the deep BEC limit, however, 
the two-body contribution to the contact arising from the 
molecular state dominates. 

This paper is structured as follows: In the next sec- 
tion, we introduce the LDA and the QMC results for the 
zero-temperature equation of state for a uniform strongly 
interacting Fermi gas. In Sec. Ill, we discuss the produc- 
tion of a strongly interacting trapped 2D Fermi gas of 
^Li atoms and how to calculate the density distribution 
and cloud size from the QMC equation of state within 
LDA. The experimental procedure for the density mea- 
surements is briefly summarized. In Sec. IV, we present 
the theoretical density distributions and sizes and com- 
pare these with the experimental measurements. In Sec. 
V. we find Tan's contact and the large- momentum static 
structure factor of a 2D trapped Fermi gas. Conclusions 
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and future perspectives are given in Sec. VI. 



II. LDA AND 2D UNIFORM FERMI GAS IN 
THE BEC-BCS CROSSOVER 

The equation of state of a strongly interacting Fermi 
gas in homogeneous space provides a convenient way to 
calculate the density distribution in a harmonic trap us- 
ing the local density approximation (LDA) [23| • The ba- 
sic idea is that for a sufficiently large number of particles 
in a slowly varying trapping potential Vext (r) we may 
treat the trapped Fermi system as a collection of many 
independent units that behave locally as a uniform Fermi 
gas. The correlation between different units, for exam- 
ple, the surface energy of each unit, is assumed to be 
negligibly small. Therefore, the local chemical potential 
of a unit at position r may be written as. 



H{v) = n [n (r)] = M - Vext (r) 



(1) 



where n (r) is the local density and ^ is the chemical po- 
tential at the trap center. At zero temperature the local 
chemical potential ji [n (r)] of the locally uniform unit de- 
pends on the local density n (r) only. Hence, given the 
local equation of state /i[n(r)] at position r, we could 
solve inversely the density n (r) = n(/i — Ve2-t(r))- The 
chemical potential at the trap center /i is set by the nor- 
malization condition / dm (r) = iV, where N is the total 
number of particles. The LDA has been shown to work 
well in a wide range of situations [23|. It is valid for ei- 
ther non-interacting or strongly interacting Fermi gases 
in different geometries from 3D to ID. 

The essential ingredient of the LDA is the local uni- 
form equation of state For a non-interacting two- 
component (spin-1/2) 2D Fermi gas at zero tempera- 
ture, the chemical potential is simply the Fermi energy 
jjL = Ep = h'^kp/{2m), where m is the mass of fermions 
and the 2D Fermi wave- vector is given hy kp = (2?™)^/^. 
Therefore, the non-interacting chemical potential is pro- 
portional to the density, /i — irfi^n/m. The mean energy 
per particle is E/N = h'^kj,/{4m) = Ep/2. 

For an interacting 2D Fermi gas in the BEC-BCS 
crossover, the energy per particle E/N has been calcu- 
lated by Bertaina and Giorgini as a function of the inter- 
action strength, by using the fixed-node diffusion Monte 
Carlo method [23. Here, we extract the chemical poten- 
tial from the QMC data of the energy per particle, since 
fi = dE/dN. A suitable parameterization of the QMC 
data for E/N is therefore needed, as we discuss in detail 
below. 

In the 2D BEC-BCS crossover, a peculiar feature of 
the contact interactions (between two fermions with un- 
like spins) is that any attraction, however small, will 
support a two-particle bound state with energy = 
-4;iV[exp(27) ma2p,], where 7 — 0.577216 is the Eu- 
ler's constant and a2D is the 2D s-wave scattering length 
[H, [Hi. This is in sharp contrast with the 3D BEC- 
BCS crossover, where a two-body bound state appears 



on only one side of the Feshbach resonance where the 3D 
scattering length is positive The scattering length 
in 2D a2D is always positive due to the existence of the 
bound state. The unitarity limit with an infinitely large 
scattering length {020 — >■ +00) is in fact trivial: it cor- 
responds simply to the non-interacting (BCS) limit. In 
the opposite (BEC) limit with infinitely small scattering 
length {a2D — >■ 0+), where the energy of the bound state 
is infinitely large, two fermions are tightly bound to form 
a composite molecule. There will be a repulsive interac- 
tion between two composite molecules, characterized by 
an effective scattering length > 0. 

The interaction strength in 2D may be expressed as 
a dimensionless interaction parameter 77 = hi{kpa2D)- 
The weakly interacting BEC and BCS limits correspond 
to 77 — > —00 and J] — > +00, respectively. The strongly 
interacting crossover regime occurs at about 77 = 0, where 
a2D ^ kp^- 

We interpolate the QMC data for the 2D equation 
of state {E/N - €3/2) /Epc in the BCS-BEC crossover 
with a smooth, continuous analytical function 7(77) that 
consists of three parts. On the BEC side in the range 
?7 < — 1/2 we use the equation of state for a gas of com- 
posite molecules with a molecular scattering length Od 
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where rud = 2m and Nd = N/2 are respectively the mass 
and number of molecules, and Ud = n/2 is their density. 
By assuming that Od/ 0,20 = Q^m ~ 0.6 [2^ , Eq. ^ turns 
into 
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On the BCS side, for 77 > 2.72 we consider a Fade-type 
approximate function 
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(4) 



where the Fade coefficients ai = 0.164106,02 = 0.702385 
and bi = 1.16411, &2 = 2.40527 are obtained by mini- 
mizing the standard deviation between the QMC data 
and the values of the function (jH) under the constraint of 
fBCs{il — ^ +00) — 1 — 1/7/, which is required in a weakly 
interacting normal Fermi liquid in 2D. In the crossover 
regime between the BEC and BCS limits we use a sixth 
order polynomial function 



/I \ ^ 

f crossover f ~ — ^7 — 2.72 j - ^ ^ C^l 



(5) 



The coefficients are selected to provide the fit best to the 
the QMC data and to ensure continuity of the equation of 
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state function and of its first and second derivatives at the 
two connection points 77 = —0.5 and 77 = 2.72. We find 
that, Co = 0.200219, cj ^ 0.154862, C2 = -0.0144822, 
C3 = 0.070831, C4 = -0.01977, C5 = 0.00891172, and 
C6 = -0.00108548. 



where lOx 



■ 10 \_ and LOz are the trapping frequencies 




Figure 1. (color online) Interpolating energy function f{r\) 
for the equation of state of a 2D Fermi gas in the BEC-BCS 
crossover. The red solid line on the BEC side is a fit corre- 
sponding to an equation of state of a gas of composite bosons, 
Eq. (|3]). The blue line on the BCS side is the Fade-type ap- 
proximation function of Eq. Q that minimizes the standard 
deviation with respect to the QMC data j23 |. The dot-dashed 
line in the crossover regime is a polynomial fit ensuring the 
continuity of the function and its first and second derivatives 
at the two connection points, see Eq. ((5}. The black circles 
show the QMC data of the equation of state. The chemi- 
cal potential derived from the interpolating energy function 
is plotted by the thin black line. 

As shown in Fig. 1, the interpolating function ](r\) 
for the equation of state provides and excellent fit to the 
QMC data. By using /i = dE/dN, we find that. 



Ef 



(6) 



The dimensionless chemical potential f^{Tf) is shown in 
Fig. 1 by the thin black hne. 



in the radial (x, y) and axial [z] directions, respectively. 
The trap aspect ratio A = ljJz/ijJi_ ^ 1. 

The basic requirement for achieving a 2D trapped 
Fermi gas may be estimated by considering first the 
zero-temperature and non-interacting limits, in which the 
properties in a 2D harmonic trap, 

Vext (r) = mwi + y^) /2 = mwipV2, (8) 

can be conveniently understood using LDA. As ^ ^ 
■nh^n/m in an ideal 2D uniform Fermi gas, we obtain that 
n(r) = mfi{r) / (Trfi^) ~ — Vext (r)]/(7rfi,^). In other 
words, we expect a Thomas-Fermi (TF) distribution. 



n(r) = fiTF 1 



PtfJ ' 



(9) 



where utf is the TF peak density and Ptf is the TF 
radius. Once the distance p is larger than the TF ra- 
dius, the density is necessarily zero. The TF peak den- 
sity and radius are related to the chemical potential by 
/i = irfi^nTF /iTT' and p, = respectively. Us- 

ing the normalization condition J dm (r) = A^, it is 
straightforward to obtain that, titf = A^^/^/(7ra5_) and 
Ptf = (4A^)^/'^a^, where aj_ = [^/(mwj,)]-^/^ is the char- 
acteristic oscillator length in the radial direction. The 
2D chemical potential or Fermi energy Ef is given by, 
p = Ef = \/NfibJ±. The characteristic Fermi temper- 
ature is Tf = EF/ks and the Fermi wave- vector fci? in 
2D harmonic traps is given by kp = \2mEF /1t?'Y^'^ = 

As the lowest excitation energy in the z direction is 
hwz, one finds that the 2D regime can be reached if 
P,Ef < hjjz and T < huz/kB- The former condition 
requires that the total number of atoms must be less 
than a 2D critical number, N2D, equal to the number 
of single particle states with energy less than the lowest 
lying state with one transverse excitation. It is straight- 
forward to show that N2D = In our experiment with 
^Li atoms, the trapping frequencies are ~ 27r x 2800 
Hz and uj± ~ 27r x 47 Hz, leading to A w 60, a± w 6.0 
pm and N2D ~ 3600. 



III. A 2D STRONGLY INTERACTING FERMI 
GAS IN HARMONIC TRAPS 

A. Achieving the 2D regime 

For an atomic Fermi gas in a harmonic trap, the 2D 
regime is achieved if the chemical potential p and tem- 
perature T are sufficiently small compared to the excita- 
tion energy in one dimension (z). We consider a spin-1/2 
Fermi gas of N ^Li atoms with equal spin-populations in 
a highly oblate harmonic trap, 

K'? (r) = l^n [iolx^ + uy + u:lz^] , (7) 



B. Achieving the strongly interacting regime 

Experimentally, the strongly interacting regime is 
reached by tuning an external magnetic field B near a 
Feshbach resonance {Bq — 834 G) for ®Li atoms, for 
which the s-wave scattering length 

a^D [B) = abg + ^^^J [l + a{B- B„)] (10) 

can be changed precisely to arbitrary values j2^. Here, 
a^g = — 1405ao with qq ~ 0.529 x 10"-^" m is the back- 
ground scattering length, AB = 300 G is the width of 
Feshbach resonance and a = 0.0004 G~^. In our highly 



4 



oblate geometry, the tight-confinement in the z-direction 
induces a bound state in the 2D x — y plane. Therefore, 
one can express the 2D scattering length in terms of the 
3D scattering length [ll.[2i|. 



a2D 



exp 



2 a^D 



(11) 



where = [^/(tow^)]^/^ and b sa 0.915. In Fig. 2, 
we plot the dimensionless interaction parameter 77 = 
ln{kpa2D) as a function of the magnetic field at the 2D 
critical number of atoms, = N2D = 3600. As indicated 
by the dotted line, at the location of the 3D Feshbach res- 
onance (B = Bo), the interaction parameter 7y w 1. 
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Figure 2. (color online) The dimensionless interaction param- 
eter ln{kFa2D) of a trapped 2D interacting Fermi gas near a 
Feshbach resonance. The vertical dotted line indicates the 
resonance position. Here, we calculate the Fermi wave-vector 
fcF at iV = N2D = 3600 by using kp = {iNy^'^aJ^. The 2D 
scattering length is calculated using Eq. pif) with ~ 770 
nm for co^ = 2tt x 2800 Hz. 



C. Theoretical density distributions 

Let us now consider the density distribution of an in- 
teracting 2D Fermi gas in the BEC-BCS crossover within 
LDA. The simple relation n{r) ~ m^{r)/ {ttIi^), useful for 
the ideal gas, is no longer applicable. We have to obtain 
numerically the local density from the local chemical po- 
tential by using the /^-function, defined in Eq. That 
is, we need to solve the following equation to find the 
density n(r) from, 

,2 



H (r) = /i - 7nujl {x^ + y^) /2, 



2 



H n (r) In y/2Tm (r)a2_D 



(12) 
(13) 



In the inversion procedure, our analytic interpolating J^- 
function appears to be very convenient. The density dis- 
tribution n(r) is calculated for an initially chosen value 



of chemical potential fi. We then adjust fi to satisfy the 
number equation J drn (r) = N. At the final stage, we 
quantify the cloud size using root mean square (rms) ra- 
dius. 



/ drn (r) 



J drn (r 
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(14) 
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Figure 3. (color online) Atomic density distributions of 
a 2D Fermi cloud at three values of interaction strengths 
77 = ln(fcFa2D). The red solid line and blue dashed line cor- 
respond to the weakly interacting BEC {rj = —3) and weakly 
interacting BCS (77 = -1-3) side. The black dotted line stands 
for a strongly interacting crossover Fermi gas (r; = 0). 

In Fig. 3, we plot the theoretical density distributions 
at three values of the interaction parameter 77, which cor- 
respond to the weakly interacting BEC and BCS sides, 
and the strongly interacting crossover regime. The den- 
sity and distance from the center of the cloud are plotted 
in units of the TF density utf and the TF radius ptf, 
respectively. One finds that with decreasing the inter- 
action parameter 77 from the BCS side to the BEC side, 
the 2D cloud becomes denser and narrower in size, as 
anticipated. On the BEC side, our analytic interpolating 
/^-function leads to the following asymptotic behavior. 
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where am = a,d/a2D ~ 0.6 is the ratio between 2D molec- 
ular and atomic scattering length. Hence, with decreas- 
ing 7/ (— > —00) the peak density increases as [—477]^/^ 
and the radius of the cloud decreases as [—477]-'^/^. On 
the other hand, on the BCS side the density distribution 
converges to the ideal Fermi gas result of Eq. (jH]), 
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(16) 
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In accord with these asymptotic density distributions, 
the cloud sizes are given by, 
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(17) 



and 
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(18) 



in the BEC and BCS Hmits, respectively. Here, 
yTp^iG — Ptf I \/3 is the rms cloud size of an ideal 2D 
Fermi gas. We have checked that these analytic results 
agree well with our numerical calculations in the appro- 
priate limits. We have also checked the sensitivity of the 
calculated density profiles to the form of the equation of 
state. The widths of the theoretical density profiles do 
not vary by more than 5% which is small on the scale of 
the width changes as the interaction strength is varied 
over the range considered here. 



D. Experimental measurements 

To measure the density distribution of a strongly in- 
teracting 2D Fermi gas, we use an experimental setup 
similar to the one used in our previous work Q . In brief, 
we start with a cloud of approximately 1 x 10^ ^Li atoms 
in the two hyperfine states \F = 1/2, mp = ±1/2) in a 
far detuned 3D optical dipole trap. The cloud is evap- 
oratively cooled to the lowest possible temperature near 
the Feshbach resonance. At this stage, the number of 
atoms is controlled by further lowering the trap depth to 
spill atoms out of the dipole trap. We then ramp on a 2D 
optical trap in 200 ms to create a highly oblate trap with 
trapping frequencies a;^ ~ 27r x 2800 Hz and uj± ~ 27r x 47 
Hz, which gives an aspect ratio of approximately 60. Fi- 
nally, we tune the interaction strength by adiabatically 
ramping the external magnetic field to 810 G, 834 G and 
992 G; where the cloud is held and imaged. The dimen- 
sionless interaction parameters hi(kFCi2D) at these fields 
are about —0.5, -1-0.6, and -|-5, respectively. 

The critical number of atoms for reaching 2D regime is 
N2D ~ 3600. Depending on the depth of the dipole trap, 
the final number of atoms in the cloud can be varied in 
the range of 500 to 5000 atoms. The final temperature 
of these small 2D and quasi-2D clouds is difficult to de- 
termine when interactions are present. However, we an- 
ticipate it to be approximately 0.1 Tp based on applying 
the same preparation procedure to clouds with a larger 
atom number. 

Before the imaging, we allow a short time of fiight 
(500 ^s). This allows us to resolve the density distri- 
bution in the tightly-confined z-direction, since this time 
scale is long compared to l/ujz ~ 57 /is. It is however 
much short compared to l/uj± « 3.4 ms, and therefore. 



the cloud distribution in the radial direction is essen- 
tially equivalent to the in situ profile. The imaging beam 
propagates roughly along the radial x-direction, which 
means that there is an automatic integration over the 
x-direction for the total density distribution. We then 
integrate these distributions over the z-direction to gen- 
erate a double- integrated column density fi{y). The rms 
cloud size is calculated from the first moment of the one- 
dimensional profile {p^} = / dyn{y)y^ / J dyh{y). The- 
oretically, we perform the same integration procedure 
in the x-direction for the 2D density distribution. This 
should lead to the same distribution as the experimen- 
tally measured double- integrated column density n{y). 



IV. COMPARISON BETWEEN THEORY AND 
EXPERIMENT 

In Fig. 4, we compare the LDA column density n{y) 
(lines) with the experimental measurements (solid cir- 
cles) at three magnetic fields. To reduce the experimental 
noise, we average the experimental density distributions 
over many images with a range of atom numbers that 
are all well below N2d- Accordingly, the theoretical lines 
correspond to the average number of atoms, while the 
standard deviation in the atom number is illustrated by 
shaded region. We observe good qualitative agreement 
between theory and experiment with no adjustable pa- 
rameters. The distributions become wider with increas- 
ing \Ti[kpa2D) however, there are notable discrepancies 
between theory and experiment, particularly in the wings 
of the clouds. This is due to a number of effects including 
the finite imaging resolution and recoil induced blurring 
during the imaging pulse. The combination of these two 
artifacts is to lower the effective resolution of the imag- 
ing system to approximately 6 /im. This alone however, 
is not enough to fully account for the observed discrep- 
ancies. The remaining differences are most likely due to 
the finite temperature of the clouds which will show up 
most in the wings of the distribution. 

In Fig. 5, we plot the rms cloud size of the 2D Fermi 
gas in units of the TF radius as a function of the in- 
teraction parameter rj = h\kpa2D- The lines and sym- 
bols show, the LDA predictions and experimental data, 
respectively. In the strongly interacting regime at the 
magnetic fields B = 810 G and 834 G, the LDA pre- 
dictions agree quantitatively well with the experimental 
data (solid red circles and green squares). At the high 
field B = 992 G, the gas is more weakly interacting and 
the experimental data (blue diamonds) lie slightly above 
the theory curve. In the inset of Fig. 5, we present the 
actual cloud size as a function of the number of atoms. 
The agreement between theory and experiment near the 
Feshbach resonance (810 G and 834 G) becomes more 
apparent. 

In all three sets of measurements, the radial cloud size 
drops below the true 2D prediction for the largest atom 
numbers. This happens when the first transverse excited 
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Figure 4. (color online) Comparison between the theoretical 
LDA predictions (lines) and experiment measurements (solid 
circles) for the column density n{y) at the magnetic fields, 810 
G (a), 834 G (b), and 992 G (c). The optical density is shown 
in arbitrary unit as a function of the actual ^-coordinate (in 
units of micrometer). 



state becomes energetically accessible and leads to a drop 
in the growth rate of the radial cloud size. In this quasi- 
2D regime, shell structure associated with resolving the 
discrete transverse states can dramatically affect the den- 
sity profiles [Hi- The atom numbers for which the radial 
width departs from the 2D theory are slightly below what 
we would expect for an ideal Fermi gas. This could be 
due to interactions or scattering resonances (27l - l3l| but 
also may arise through thermal excitations. 




Figure 5. (color online) The rms cloud size in units of the TF 
radius as a function of the dimensionless interaction strength. 
The black line shows the theoretical LDA prediction. The 
red circles, green squares, and blue diamonds represent the 
experimental data at the magnetic fields, 810 G, 834 G, and 
992 G, respectively. The number of atoms at a given magnetic 
field can be varied, giving rise to slightly different interaction 
parameter. In the inset, the actual rms cloud sizes are shown 
as a function of the decimal logarithm of the atom number. 
The blue solid line (the highest) is the LDA prediction for 992 
G; the green dashed line (medium) is for 834 G; and the red 
dot-dashed line (the lowest) is for 810 G. The vertical dotted 
line in the inset indicates the position of the critical atom 
number. 



V. CONTACT PARAMETER AND THE STATIC 
STRUCTURE FACTOR 

A. Tan's contact parameter 

For a strongly interacting Fermi gas with contact in- 
teractions, the asymptotic behavior of various physi- 
cal quantities in the limit of short-distance or large- 
momentum is governed by a single parameter, called the 
contact, which measures the density of fermionic pairs 
within a short distance. This was first discussed by Tan 
in 2008, when he derived a set of exact universal rela- 
tions for strongly interacting Fermi gases [23j . Being an 
important many-body parameter. Tan's contact is also 
related to the thermodynamics via the adiabatic relation 
(23| . The contact in a homogeneous 2D Fermi gas was 
calculated by Bertaina and Giorgini Here we use 

these results to find the contact and static structure fac- 
tor in a trapped system. In 2D, the adiabatic relation 
takes the form [331, 



27rm dE 
d In a2D ' 



(19) 



where the derivative is taken at constant entropy. The 
calculation of Tan's contact parameter for a 3D strongly 
interacting Fermi gas was performed recently (33| , by us- 
ing the similar adiabatic relation. 
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The zero-temperature contact of a homogeneous 
2D Fermi gas can be calculated by substituting our 
interpolated energy per particle E/N = 6^/2 + 
EFGf{^^[kFa2D\) mto Eq. We find that I = 

l2b +Imb, where 



1-2b 

Nkl 



27rm d[ei3/2] 



Stt 1 



(20) 



is the contribution from the two-body bound state and 
Trnb Sttto dlEpcf iv)] 



Nkl 



d\na2D 



2 dr] 



(21) 



is the contribution from the many-body correlations, re- 
spectively. The two-body contact l2b increases mono- 
tonically from the BCS to the BEC limit. In contrast, 
the many-body contact Xmb should exhibit a maximum 
at the crossover regime, according to the behavior of the 
energy /-function (see Fig. 1). As shown in Fig. 6 by a 
thin line, the maximum of the many-body contact occurs 
at 77 ~ 0.8, which is roughly the position of the Feshbach 
resonance (see also, ref. |23|). 




Figure 6. (color online) Theoretical contact parameter of a 
trapped 2D Fermi gas in the BEC-BCS crossover. The con- 
tact is shown in units of Nk%. The blue dot-dashed line and 
red dashed line show respectively the contribution from the 
two-body bound state and the many-body part (see text), 
while the thick black line gives the total contribution. For 
comparison, we show by a thin line for the many-body part 
of the contact of a homogeneous 2D Fermi gas. 



However, the many-body part may be significantly af- 
fected. In Fig. 6, we present the result for the contact of 
a trapped 2D Fermi gas at the BEC-BCS crossover, with 
the many-body contact shown by a red dashed line. It 
has roughly the same shape as the many-body contact of 
a homogeneous gas, with a peak appearing at ry ^ 0.8. 
However, the peak is about half as high due to the aver- 
age over the density distribution. 



B. Spin-antiparallel static structure factor 

Tan's contact for a 3D strongly interacting Fermi gas 
has been measured in a number of ways. One appealing 
method is to measure the spin-antiparallel static struc- 
ture factor by using Bragg spectroscopy at large momen- 
tum [s^l, which has a 1/q tail with a prefactor given by 
Tan's contact |35| . In 2D, we may make a similar pre- 
diction. It has been shown that the 2D pair correlation 
function n'^) (r) cx Iln^(p/a2r,) HI- The ln^(p/a2D) de- 
pendence can be qualitatively understood from the two- 
body relative wave- function -iprei (r) ~ ln(r/a2_D), since 
n'^^ (r) oc \iprei (r)l^. The spin-antiparallel static struc- 
ture factor is simply the Fourier transform of the pair 
correlation function. Thus, we find that. 



Sn (g > ^f) 



I [7 + ln(ga2_D)] 



N 



9 9 



Nkl 



[7 + In {kFa2D) + Ing] 



(23) 



,(24) 



where q = q/kp and 7 = 0.577216. Compared with 
the 3D case, the spin-antiparallel static structure factor 
in 2D decays faster with increasing momentum q (9^^ 
compared to q"^). 

In Fig. 7 we plot the theoretical prediction for the spin- 
antiparallel static structure factor of a trapped 2D Fermi 
gas in the BEC-BCS crossover for momenta q = Skp and 
q = bkp- We split the structure factor into the two-body 
and many-body parts, in accord with the previous clas- 
sification for the contact. Near the Feshbach resonance, 
the static structure factor a.t q ^ Skp is about 0.02, whose 
magnitude is accessible within current experimental res- 
olution o. 



VI. CONCLUSIONS 



For a trapped interacting Fermi gas, we may calculate 
the contact by using the LDA density distribution. 



dr 



I(r) 
AV 



2n 





I 


j dr 




[Nkll 



(r)n2(r), (22) 



where we have summed the local contact density, 
I{r)/AV = [I/{Nkl)] (r) x n(r)fc|,(r), over the whole 
space. It is easy to see that, the two-body contact is 
not affected by the density average, so that It.26 = 2^26- 



To conclude, we have predicted theoretically and mea- 
sured experimentally the density distribution and cloud 
size of a low-temperature two-dimensional harmonically 
trapped Fermi gas in the BEC-BCS crossover. The the- 
oretical calculations have been carried out within a lo- 
cal density approximation, based on the ab-initio zero- 
temperature equation of state obtained from the fixed- 
node diffusion Monte Carlo simulations . The experi- 
mental measurements were performed using a single two- 
dimensional Fermi cloud of ^Li atoms near a Feshbach 
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resonance. We have found good qualitative agreement 
between theory and experiment. 

We have also calculated an important many-body pa- 
rameter, Tan's contact, and have proposed that it can 
be straightforwardly measured using Bragg spectroscopy 
for the spin-antiparallel static structure factor as in three 
dimensions. In future studies, it will be interesting to 
study both experimentally and theoretically the density 
distributions at finite temperatures, which may elucidate 
the fermionic Berezinskii-Kosterlitz-Thouless transition 
in two dimensions. 



Figure 7. (color online) The LDA prediction for the spin- 
antiparallel static structure factor of a trapped 2D Fermi gas 
in the BEC-BCS crossover, under our experimental conditions 
of ~ 27r X 2800 Hz and u]^c^2n x 47 Hz for iV = N2D = 
3600 ^Li atoms. We show the many-body contribution by the 
dashed lines. The thick and thin lines give respectively the 
predictions at the transferred momentum g = 5A;f and 3A;f. 
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